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Abstract 

Non- vacuum exact gravitational waves invariant for a non Abelian two- 
dimensional Lie algebra generated by two Killing fields whose commutator 
is of light type, are described. The polarization of these waves, already 
known from previous works, is related to their natural sources consisting 
of cosmic strings and 7-ray bursts. Non vacuum exact gravitational waves 
admitting only one Killing field of light type are also discussed. 

PACS numbers: 04.20.-q, 04.20. Gz, 04.20. Jb . 



Introduction 

Because of experimental efforts devoted to the detection of gravitational waves, 
there is a growing theoretical interest in describing the emission of detectable 
gravitational waves from astrophysical sources. All experimental devices, laser 
interferometers or resonant antennas, need very precise theoretical predictions, 
and usually they are calibrated according to results of linearized theory. In fact, 
it is commonly believed that, since the sources are very far from the earth, the 
wave's amplitude is so small at the earth surface that the use of the linearized 
theory is lawful. However, in some cases this assumption is misleading. Indeed, 
as shown by Christodoulou [3] the nonlinearity of the Einstein equations mani- 
fests itself in a permanent displacement of the test masses of a laser interferom- 
eter which, in many physically interesting cases (e.g. for a binary coalescence), 
cannot be neglected. Moreover, the nonlinearity of the Einstein theory is one of 
its most characteristic feature, and it is likely that some crucial properties of the 
gravitational field are encoded in the nonlinear terms. Thus, a detailed analysis 
of exact gravitational waves, besides to help the calibration of the experimental 
devices, could disclose new phenomena hidden by the approximations. 



1 



In recent papers exact solutions [14] of vacuum Einstein field equations, 
representing spin-1 gravitational waves [4, 5], have been discussed. They are 
described by the following metrics: 

g = 2f (dx 1 dx 1 + dx 2 dx 2 ) + n [(w(x\x 2 ) - 2x 4 ) dx 3 dx 3 + 2dx 3 dx 4 ] 

(1) 

where, [i = Z)$ + B; D, B e 1Z, $ is a non constant harmonic function of 
a; 1 and x 2 , f = (V$) 2 a/|^|//z, and w (a; 1 , a; 2 ) is a non constant 1 solution of the 
/i- deformed Laplace equation: 

Aw + (d x i In |/z|) d x iw + (d x 2 In |^|) d x 2 W = 0, (2) 

which, when /j, is a constant function, reduces to the Laplace equation. 

These metrics are intrinsically characterized [14] by the following properties: 

• They are invariant for the non Abelian two-dimensional Lie algebra gen- 
erated by the Killing vector fields X and Y such that: 

[X, Y]=Y, 



with Y of light-type, i.e., g (Y, Y) = 0. 

• The distribution T> , orthogonal to the distribution V generated by X 
and Y, is integrable and transversal to T>. 

In the following, an integral (two-dimensional) submanifold S of T> will be 
called a Killing leaf and an integral (two-dimensional) submanifold W of V 1 - 
will be called a orthogonal leaf. 

In [14] it is shown that these Einstein manifolds have a have a natural fibered 
structure with the Killing leaves as fibers and a basis diffeomorphic to the or- 
thogonal laves. 

In the coordinates {x^), which are adapted [14] to the Killing fields, the 
Killing vector fields X and Y read 

In the particular case / = 1/2 and \i = 1, the above metrics are locally 
diffeomorphic to a subclass of the vacuum Peres solutions [8] corresponding to 
a special choice of the harmonic function parameterising that metrics. 

This is easily seen by introducing, for it ^ 0, new coordinates (x,y,u,v) 
defined by 

x 1 = x, x 2 = y, x 3 = In \u\ , x 4 — uv, 



1 Whcn w is constant, metrics (1) are static. 
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in which metrics (1) read 



g = 2f(x, y) \dx 2 + dy 2 ] + p,(x, y) 2dudv H — ^du 2 . 



(3) 



A crucial problem in the theory of the gravitational waves is to characterise 
realistic sources able to generate waves so strong to be detected by the exper- 
imental devices. Thus, the aim of this paper is to find physically interesting 
sources for the above metrics and to relate the physical properties of such grav- 
itational waves (in particular, the asymptotic flatness, the wave-like character 
and the polarization) to the sources. The plan of the paper is the following: in 
the first section, we will solve the non- vacuum Einstein equations for metrics (1) 
with a dust source describing a null electromagnetic wave. In the second section 
we will discuss the role of the asymptotic flatness and the related constraints 
on the metric components and on the source. In this section we will also dis- 
cuss cosmic strings as a possible source. Then we will briefly explain why the 
standard textbook analysis of gravitational waves does not apply to this case. 
In the third section we will show that such gravitational fields indeed represent 
gravitational waves by using two well known covariant criteria: the Zcl'manov's 
and the Pirani's criteria. In the fourth section we will see that such waves are 
spin-1 objects. In the last section we will give a hopefully vague discussion on 
the possibility to detect such non standard gravitational waves. 

1 The dust and cosmic string sources 

The simplest source for metrics (1) is dust with density p and velocity £/ M and, 
then, characterized by an energy-momentum tensor T Mi , = pU^U^. When /7 M 
is a light-like vector field, this kind of energy-momentum tensor can describe 
the energy and momentum of null electromagnetic waves, i.e., electromagnetic 
fields whose scalars, F^ v F^ and e a/3f11 ' F a pF^, both vanish. In this way it could 
be possible to describe the gravitational effects (in particular the emission of 
gravitational waves) of a very interesting astrophysical phenomenon, the 7— ray 
bursts (GRBs): emission of ultra-high energetic 7-rays, whose origin is still to 
be fully understood. 

Moreover, with this source one can also preserve the symmetries [14, 4, 5] of 
the vacuum solution. 

Being the time coordinate in the Killing leaves, the dust cannot move orthog- 
onally to them and it will be chosen to move parallel to the light-like Killing 
field Y, i.e., with velocity = S^ 4 . 

Thus, the metric will be taken in the form 



<) 




+M(x 1 ,x 2 ) \w(x\x 2 ,x 4 ) (dx 3 ) 2 + 2dx 3 dx 4 



(4) 
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with F, M and W to be determined by the non vacuum Einstein equations 
with the following energy-momentum tensor 

V = MV„3<W (5) 

Looking at the explicit expression of the Einstein tensor it turns out 
that 

G 33 = WG 34 + -5- (MAW + VM • VW) , 
It 

where 

9 9 \ „„ „ T „ dikfcW 9M9W 
VM- Vlf = — -— - + 



9a; 1 ' 9x 2 / ' 9a; 1 9a; 1 9a; 2 9a; 2 

Then, since = for u,, z/ 7^ 3, we can fulfil the non vacuum field equations, 
while keeping the vacuum Killing fields, if and only if F = / and M = fi, 
i.e. F and M coincide with their vacuum expression [14], and W has the form 
W (x 1 , a; 2 , a; 4 ) = w (x 1 , x 2 ) — 2x 4 , w being a solution of the [i- deformed Poisson 
equation 

fxAw + V/i • Vw = 2//x 2 p, (6) 

where, to save writing, beyond c = 1, it has been taken 8ttG = 1, G being the 
Newton gravitational constant. 

The comparison of Eqs. (2) and (6) shows that the introduction of the 
matter source (5) induces only a simple and natural modification in the field 
equations. The reason behind this simplicity is that the above source, besides 
being physically interesting, does not break the geometrical structures which 
allows [14] to classify completely Ricci-flat metrics of the form (1). In particular, 
in this case it is still true that the 2-dimensional distributions orthogonal to the 
surface spanned by the Killing fields is integrable since the energy-momentum 
tensor has no components orthogonal to the Killing leaves. Thus, the bundle 
structure [14] of the vacuum space-time is not destroyed by the source. 



2 Asymptotic flatness 

From the physical point of view, it is important to characterise, among the 
metrics (1), those which are spatially asymptotically flat. The vacuum case, 
where (for fi = 1) x 1 and x 2 are harmonic coordinates, suggests to call (spatially) 
asymptotically flat a metric approaching, for (x 1 ) + (a; 2 ) — > 00, the Minkowski 
metric. 

This intuitive definition of asymptotic flatness allows to obtain qualitative 
results by using the standard theory of partial differential equations. 

In terms of the functions /, \i and w, the asymptotic flatness condition reads: 

(x 1 ) 2 + (a; 2 ) 2 -> 00^ 

/ — > const, yU — > const, w — > cix 1 + c 2 x 2 + C3, 
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where c\ , ci and C3 are arbitrary constants and the behaviour of w can be easily 
recognized by looking at the Riemann tensor. 

It is easy to see that, in order for metrics (1) to be spatially asymptotically 
flat, necessarily p must be constant. Indeed, being / and p related by 

/ = ( V $) 2 , p = A$ + B, A$ = 0, 

M 

the condition p — ► const implies A = or $ — > const. However, the case 
$ — ► const has to be excluded because of the nondegeneracy of the metric. 

For p = 1, the equation for to reduces to the two-dimensional Poisson equa- 
tion 

Aw = p. 

It is well known that, if p goes to zero fast enough, it is possible to find non 
trivial everywhere regular solutions w tending to a constant value 2 . 
The function / satisfies the equation 

/A/ - (V/) 2 = 0, 

and this implies that the function 

V = ln|/| 

is harmonic. Therefore, to have / asymptotically constant, the function ip must 
be asymptotically constant. But this is impossible unless ip = const. Thus, in 
order to have everywhere regular spatially asymptotically flat solutions, / and 
p must be constant functions and the fluid density p must tend to zero fast 
enough. 

However, if we admit 5-like singularities in the (a; 1 ,^ 2 ) plane (i.e., string- 
like singularities, by taking into account the third spatial dimension), spatially 
asymptotically flat vacuum solutions with / ^ const and w ^ const can exist. 
It is worth noting that, in this limiting case in which ^(a; 1 ,^ 2 ) — > S (a; 1 ,^ 2 ), 
the energy-momentum tensor becomes the one usually employed to describe the 
gravitational effects of topological defects known as cosmic strings [15]. This 
kind of extended objects are predicted in some particles physics cosmological 
models with phase transitions. Moreover, cosmic strings could have an impor- 
tant role in the description of two very interesting, but not yet fully understood, 
astrophysical phenomena: the GRBs and the ultra high energy cosmic rays, i.e. 
cosmic rays with energy E > 10 11 GeV, (for a review, see, for example, [13]). 
It is important to stress here that cosmic strings are not singular objects and 
should be described by a regular energy-momentum tensor, but the above widely 
used scheme is useful for practical computation. 

2 In the vacuum case this is not possible, unless w = const, because the Laplace equation 
has not solution of this kind, as can be easily checked, for example, by using the maximum 
principle. 
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Thus, from the phcnomcnological point of view, it is very intriguing that 
these kind of (wave-like, as we will show later on) gravitational fields are nat- 
urally coupled to cosmic strings. Eventually, unlike the electromagnetic (or 
the standard gravitational) waves that can be everywhere regular even without 
sources, these gravitational fields, to be globally nonsingular, have to be cou- 
pled with matter sources. To be more precise, the reason why in the standard 
textbook analysis of the polarization spin— 1 gravitational waves always appear 
as pure gauges is that one usually considers only square-integrable perturba- 
tions (in fact the Fourier analysis is needed to kill such spin— 1 components) to 
the Minkowski metric. However it is easy to verify that, if one considers non 
square-integrable perturbations, the spin— 1 components cannot be gauged away 
anymore [4]. In the case we are here considering, the physical components of the 
gravitational field arc not square-integrable and represent (as we will see later 
on) spin— 1 gravitational waves. 



3 Wave-character of the field 

From now on, to have a clearer physical interpretation, only the case / = 1/2 
and /i — 1 will be considered. In fact, in this case the metric (1), besides being 
an exact solutions of the non vacuum Einstein equations, is a solution of the 
linearized Einstein equations on flat background too. 

As it has been said, metrics (1) and (3) are locally diffcomorphic. Moreover, 
for f = 1/2 and /i = 1, the metric (3) is a particular case of the Peres metric 
[8]: 

ds 2 = dx 2 + dy 2 + 2dudv + h(x, y, u)du 2 (7) 

where h is an arbitrary function of u. If h xx + h yy — 0, the metric (7) is a solution 
of vacuum Einstein field equation, otherwise, since one still has g^R^v = and 
g^ v R^aR^v = 0, the source can be interpreted as a null electromagnetic field 
[7, 9]. 

With the chosen energy-momentum tensor, the non-vacuum Einstein equa- 
tions for the metric (7) reduce to 

A/ l = A4 = 4. (8) 

Now, let us briefly analyze the non vacuum Peres metric. 

Firstly, it is trivial to show that (x,y,u,v) are harmonic coordinates. Fur- 
thermore, because of the Laplacian operator in Eq. (8), the vacuum (i.e. with 
p = 0) Peres solution also can be asymptotically flat only if we admit S— like 
singularities in the (x, y) plane. In other words, the vacuum Peres solution is 
naturally coupled to cosmic strings too. Eventually, it is easy to see that, in 
general, this metric has only one 3 (null) Killing field: V = d/dv and that the 

3 Of course, if h = w/u 2 , then there exist an additional Killing field. 
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non vanishing independent components of the Ricmann tensor are: 

The wave character and the polarization of these gravitational fields can be 
analysed in many ways. For example, we could use the Zel'manov criterion 
[16] to show that these are gravitational waves and the Landau-Lifshitz pseudo- 
tensor to find the propagation direction of the waves [4, 5]. However, the alge- 
braic Pirani criterion is easier to handle since it determines the wave character 
of the solutions and the propagation direction both at once. Moreover, it has 
been shown that, in the vacuum case, the two methods agree [5]. To use this 
criterion the Weyl scalars must be evaluated according to the Petrov-Penrose 
classification [10, 11]. 

To perform the Petrov-Penrose classification, one has to choose a tetrad basis 
with two real null vector fields and two real spacelike (or two complex null) 
vector fields. Then, according to the Pirani's criterion, if the metric belongs to 
type N of the Petrov classification, it is a gravitational wave propagating along 
one of the two real null vector fields. However, the Peres coordinates are not of 
the above type and it is necessary to introduce new coordinates adapted to the 
Petrov-Penrose classification. 

By performing the following transformation 

x i ► x , y i — ► y, mhii, v v + ip (x,y,u) , 

whose Jacobian is equal to one, the metric (7) with the choice h — <p iU becomes: 

g = dx 2 + dy 2 + 2dudv + 2(tp tX dx + ip^ y dy)du (9) 

and has, for generic ip, the only Killing field Y = d/dv. 

Since d u and d v are null real vector fields and d x and d y are spacelike real 
vector fields, the above set of coordinates is the right one to apply for the Pirani's 
criterion. 

Since the only nonvanishing components of the Riemann tensor, correspond- 
ing to the metric (9), are 

Riuju — dijd u (p, i,j — x, y, 

this gravitational fields belong to Petrov type N ([2, 16]). Then, according 
to Pirani's criterion, the metric (9) does indeed represent a gravitational wave 
propagating along the null vector field d u . It is worth noting here that, in 
order to have a spatially asymptotically flat gravitational field, d u ip has to be a 
solution of the two-dimensional Poisson equation in approaching a constant at 
infinity; therefore, d u ip, ip yX and ip tV are not square-integrable functions. 
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4 Polarization 



Once the propagation direction has been determined, to compute the polar- 
ization we only need to look at the transformation properties of the physical 
components of the metric (9) under a rotation in the (x, y) plane orthogonal 
to the propagation direction. By "physical components" we mean the true de- 
grees of freedom of the metric. Since the metric (9) has only two independent 
components, we have simply to look at the transformation properties of these 
components. Moreover, Bel's superenergy tensor singles out the same degrees of 
freedom [5]. It is worth noting here that the metric (9), besides being an exact 
solution of the Einstein equations, is also a solution of the linearized Einstein 
equations on flat background: 

g = dx 2 + dy 2 + 2dudv + 2(ip tX dx + (p^ydy)du = rj a b + h a b 

hab = S au Sbif.i 

where r\ a b is the flat metric and h a b can be considered as the " perturbation" . 
Then, the spin may be discussed using the standard textbooks analysis. How- 
ever, before working out the transformation properties of the physical compo- 
nents of the metric, we should verify that the system of coordinate is harmonic 
(otherwise, there would be no meaningful notion of polarization). In fact, by 
construction, the system of coordinates of the metric (9) is harmonic, so that 
the standard procedure applies. Thus the analysis of the polarization of this 
gravitational waves runs in exact the same way as in [5]. The result is the same: 
these are spin-1 gravitational waves because, roughly speaking, the (non square- 
integrable) physical components of these waves h ux — ip, x and h uy — <p iy have 
only one index in the (x,y) plane orthogonal to the propagation direction d u . 
It is worth noting that these results also hold when the function / in Eq. (1) 
is not constant. The reason behind this non-standard result could be that the 
spacetime topology is non trivial due to the singularities of the function d u <p (as 
can be easily recognized by considering the case in which d u (p is a solution of 
the Poisson equation with S— like sources, or, in other words, when we consider 
the physical effects of cosmic strings). The gauge transformations are able to 
kill spin-1 degrees of freedom only locally, but in this case the singularities seem 
to be a global obstruction. 

The possibility of spin— 1 gravitational waves has also been guessed, although 
in a rather different context, in [1] and [6]. 

5 Possible experimental tests 

The following is a rather vague discussion on the possibility to detect such non 
standard gravitational waves. The vagueness of the below comments is mainly 
related to the lacking of experimental data on gravitational waves. Nevertheless, 
we believe that it is interesting to outline in a very general (vague) way how 
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one could recognize the presence of spin-1 gravitational waves in a would be 
gravitational wave's bursts detected by a laser interferometer. The effects of 
these waves on neighbouring test particles can be addressed by using the Jacobi 
geodesic deviation equation which for the metric (9) reads: 

= -g^djdh^, (10) 

where Z 11 is a spacelike vector, with components only in the (x, y) plane, repre- 
senting the displacement of two close geodetics and r is the proper time along 
the close geodesies. 

The above equation shows that one can have either attraction or repulsion 
according to the choice of <p, this <p(x, y, u) being constrained, outside the matter 
source, only by the condition to be a harmonic function of x and y. For example, 
the choice ip = p(x,y)a(u), with p a harmonic function of x and y and a a 
decreasing function of u, will give repulsion. This is not surprising because it is 
known from QFT that spin-odd bosons generate repulsion between particles of 
the same charge, the charge, in this case, corresponding to the mass. 

Now, the question is: "Can we hope to detect such spin-1 waves?". From 
the theoretical point of view, it would be important to have an estimate of the 
flux of spin-1 gravitational waves. However, it is rather difficult to obtain even 
a rough estimate of the energy-density of all possible sources of such waves. For 
this reason an approximate expression for p is not available yet, and, without 
p, it is not easy to estimate the expected flux. From the experimental point 
of view, as it is well known, there are delicate technical problems to detect 
gravitational waves from astrophysical sources because of the weakness of the 
wave amplitudes with respect to experimental noises. Nevertheless, once it will 
be possible to "observe" gravitational waves, it could be possible to distinguish 
the spin-2 and the spin-1 components. Indeed, suppose to put the detector (such 
as a laser interferometer or a resonant antenna) in the (x, y) plane orthogonal to 
the propagation direction of the waves. By definition, if a spin-2 gravitational 
wave's burst is observed, then the experimental data have to be invariant under 
a rotation of 7r in the (x, y) plane. This means that the observable effects of 
the gravitational wave's burst do not change if the detector is rotate of an angle 
equal to tt. However, if the burst contains also a spin-1 component, then the 
above conclusion is not true anymore. In fact, the observable effects of the spin- 
1 component are invariant under a rotation of 27r in the (x, y) plane. Therefore, 
a breaking of the invariance of the experimental data under a rotation of 7r in 
the (x,y) plane could be interpreted as the effect of spin-1 gravitational waves 
and the amount of such a breaking could be used to get an estimate of the 
flux of such waves. Moreover, since cosmic strings are naturally coupled to 
spin-1 gravitational waves, if the device is able to detect mainly cosmological 
gravitational waves, it could be possible, by looking at the experimental data, 

4 The choice p = In \f x 2 + j/ 2 gives an asymptotic flat solution. 
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to put further constraints on the presence of cosmic strings. In fact, even a 
small spin-1 component could be the signature of a non zero density of cosmic 
strings in the early cosmology. 

References 

[1] D. V. Ahluwalia, N. Dadhich, M. Kirchbach, Int. J. Mod. Phys. D 
(2002)1621. 

[2] S. Chandrasekar, The mathematical theory of black holes, (Clarendon Press, 
Oxford, 1983). 

[3] D. Christodoulou, Phys. Rev. Lett. 67, 1486 (1991). 

[4] F. Canfora, G. Vilasi, P. Vitalc, Phys. Lett. B, 545 (2002)373-378. 

[5] F. Canfora, G. Vilasi, P. Vitale, Spin-1 gravitational waves, gr-qc/0212024; 
to appear in Int. J. Mod. Phys. B. 

[6] M. Kirchbach, D. V. Ahluwalia, Phys. Lett. B, 529 (2002) 124-131. 

[7] C. W. Misner, J. A. Weeler, Ann. Phys. (N.Y.) 2, 525 (1957) 

[8] A. Peres, Phys. Rev. Lett. 3, 571 (1959) 

[9] A. Peres, Phys. Rev. 118, 1105 (1960) 

[10] A. Z. Petrov, Einstein spaces, (Pergamon Press, New York 1969) 

[11] R. Penrose, Ann, of Phys. 10:171, (1960) 

[12] N. Rosen, Bull. Res. Coun. Isr. 3, 328 (1954). 

[13] M. Sakellariadou, hep-ph/0212365 (2002). 

[14] G. Sparano, G. Vilasi and A.M.Vinogradov, Phys. Lett. B 513, 142 (2001); 
Diff. Geom. Appl. 16, 95 (2002); 17, 1 (2002). 

[15] A. Vilenkin, E. P. S. Shellard, Cosmic Strings and Other Topological De- 
fects, (Cambridge University Press, 2000). 

[16] V. D. Zakharov, Gravitational waves in Einstein's theory, (Halsted Press, 
N.Y.1973) 



10 



